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1  INTRODUCTION 

Moat  materials  behave  In  a  brittle  manner  particularly  If  operating 
temperatures  are  low  and  rates  of  loading  are  high.  However,  there  are  many 
materials  which  are  brittle  even  at  ambient  temperatures  and  slow  rates  of 
loading.  Typical  examples  are;  most  engineering  ceramics,  tungsten  below  the 
ductile-brittle  transition  temperature,  graphites,  carbons,  concrete,  cement 
and  compressed  powder  compounds.  The  characterisation  of  these  materials, 
the  manufacture  of  components  and  their  response  to  load  require  special 
treatment. 

The  measurement  of  the  strength  of  a  material  Is  usually  made  for  one  of  two 
purposes;  either  to  compare  the  relative  merits  of  two  materials  or  to 
provide  design  data  for  structural  components.  If  the  material  Is  ductile 
the  stress  concentrations  which  may  arise  In  the  test  pieces  are 
redistributed  through  local  plastic  flow.  A  brittle  material  Is  far  less 
forgiving  In  Its  response  to  load  and  if  a  flaw  becomes  critical  It  will 
usually  propagate  and  result  In  catastrophic  failure.  Mechanical  testing  of 
brittle  materials  requires  accurate  specimen  alignment  and  Its  neglect  can 
cause  premature  failure  leading  to  spurious  test  results.  This  Is 
particularly  true  for  the  simple  tensile  test  and  even  though  It  would 
simplify  the  numerical  analysis  It  Is  rarely  chosen  for  brittle  material 
characterisation.  More  suitable  test  configurations  In  which  alignment  Is 
not  so  critical  are  the  flexure  of  beams  and  the  annular  bending  or  diametric 
compression  of  discs. 

For  brittle  materials  It  Is  common  to  observe  a  large  scatter  In  the  fracture 
data  of  nominally  Identical  teat  pieces.  This  Is  due  to  the  fact  that  volume 
and  surface  defects  are  randomly  distributed  throughout  the  material  and  vary 
In  occurrence,  orientation  and  severity.  Compared  with  that  of  most  other 
structural  materials,  the  variability  of  strength  of  brittle  materials  Is 
high.  This  is  reflected  in  the  analysis  of  strength  test  data  where  large 
coefficients  of  variation  are  observed  -  typically  15Z-20Z,  and  values  as 
high  as  30Z  are  not  uncommon.  These  figures  Imply  that  a  deterministic 
strength  and  an  associated  factor  of  safety  are  no  longer  sensible  design 
parameters.  Instead  a  statistical  approach  Is  required  to  account  for  the 
material  variability. 

The  distribution  function  generally  adopted  In  brittle  material  design  Is  the 
one  due  to  Welbull.  This  distribution  Is  derived  using  a  weak  link 
hypothesis  and  Its  use  allows  the  component  failure  probability  to  be 
expressed  In  terms  of  the  failure  characteristics  of  an  elemental  volume  or 
surface.  The  Welbull  equation  can  also  model  fracture  data  with  skewed 
distributions  and  If  the  three  parameter  equation  Is  adopted  an  allowance  can 
be  made  for  a  non-zero  threshold  strength,  le  a  strength  below  which  the 
probability  of  failure  Is  zero. 

The  weak  link  hypothesis  leads  to  a  sice  effect  relationship  In  which  large 
components  are  predicted  to  be  weaker  than  smaller  ones  since  It  Is  more 
likely  that  a  critical  flaw  will  be  encountered  as  the  volume  or  surface  area 
increases.  For  materials  with  low  variability  (Is  most  metals),  the  sice 
effect  la  negligible  and  strength  data  are  not  quoted  with  respect  to 
component  sice.  However,  for  brittle  materials  where  the  variability  Is 
significant,  strength  must  be  related  to  either  a  unit  volume  or  a  unit  area 
depending  on  Aether  volume  or  surface  defects  are  considered  significant. 

The  characterisation  of  brittle  materials  requires  careful  experimentation 
and  la  normally  achieved  by  analysing  the  fracture  data  obtained  from 
flexural  testa  on  alender  beams.  The  same  geometry  can  be  used  to  study  both 
volume  and  surface  defects  providing  the  appropriate  analysis  Is  applied  In 
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each  case.  Once  the  characteristic  failure  properties  of  a  material  have 
been  evaluated,  the  failure  probability  of  a  component  subject  to  a 
aultiaxlal  stress  state  can  be  examined.  The  development  of  probabilistic 
failure  theories  which  allow  this  have  been  the  subject  of  considerable 
attention  during  recent  years. 

In  section  2  of  this  report  the  three  parameter  Weibull  equation  Is 
Introduced  and  Its  applications  to  brittle  naterlal  design  are  described. 
Section  3  presents  a  new  method  for  the  evaluation  of  the  Helbull  parameters 
based  upon  a  non-linear  least  squares  analysis.  The  generalisation  of  the 
simple  Weibull  equation  for  non-uniform  stress  fields  and  Its  subsequent 

application  to  the  analysis  of  volume  and  surface  defects  are  given  in 

Section  4.  Section  5  gives  details  of  a  multlaxlal  stress  theory  for 
predicting  the  failure  probability  of  certain  types  of  brittle  materials. 

In  another  report  by  the  authors  (Ref  1),  various  experimental  results  are 
analysed  and  the  statistical  analysis  developed  In  this  work  Is  assessed. 

2  WEIBULL  DISTRIBUTION 

From  a  given  set  of  fracture  data  the  frequency  distribution  of  strength  can 
be  constructed  by  recording  the  number  of  failures  at  the  various  stress 

levels.  For  many  materials  the  frequency  distribution  is  strongly  skewed  and 

the  theoretical  distribution  chosen  for  any  subsequent  analysis  must  be 
capable  of  describing  this.  It  Is  also  observed  that  there  is  a  threshold 
strength  below  which  no  component  will  fall  and  this  Implies  the  existence  of 
a  lower  positive  bound,  below  which  the  probability  of  failure  Is  zero. 

These  two  features  Immediately  rule  out  the  normal  distribution  which  Is 
symmetric  and  Is  defined  over  a  range  extending  from  minus  to  plus  infinity 
(Ref  2).  Using  a  normal  distribution  would  result  In  a  finite  probability  of 
failure  value  for  a  component  under  zero  load. 


A  distribution  which  has  been  used  extensively  to  characterise  brittle 
materials  Is  the  one  due  to  Weibull  (Ref  3).  The  frequency  distribution  for 
the  three  parameter  equation  In  terms  of  an  applied  stress  o  Is 


F  -  0  OiOj, 


and  the  corresponding  cumulative 
equation  r  -i 


P  -  1 


P  -  0 


probability  distribution  is  given  by  the 
oxjj., 


OtOj. 


(2) 


The  Weibull  parameters  m,  oQ  and  Oj,  are  obtained  by  fitting  the  above 
distribution  to  the  experimental  data.  In  Equation  2  the  threshold  stress 
otf  defines  the  stress  below  which  the  probability  of  fal.'ure  Is  sero. 
Frequency  and  cumulative  distribution  curves  with  m  ■  3  and  cQ  ■  1  for  ■ 
0.0,  0.3,  1.0  are  shown  In  Figure  1.  Changing  the  threshold  value  while 
keeping  the  other  two  parameters  fixed  results  In  s  simple  translation  of  the 
distribution  curves.  The  physical  meaning  of  the  Weibull  modulus  m  can  be 
seen  from  Figure  2  where  the  frequency  and  cumulative  probability 
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distributions  are  plotted  for  oT  ■  0,  Oq  •  1  and  a  -  2,  3,  4.  As  a  Increases 
the  spread  of  the  frequency  distribution  decreases,  becoming  a  single  spike 
as  a  tends  to  «•.  This  corresponds  to  a  step  change  In  the  cuaulatlve 
distribution.  Syaaetrlc  distributions,  eg  the  noraal  distribution,  are  said 
to  have  zero  skewness  and  this  occurs  In  the  Welbull  distribution  when  a  “ 
3.602.  It  Is  also  Interesting  to  note  that  all  the  probability  curves  In 
Figure  2  pass  through  the  saae  point  given  by  the  coordinates 


0  "  °T  +  °o  * 
P  -  1  -  1/e  • 


(3) 


The  noraallslng  parameter  a  affects  the  'scale'  of  the  distribution.  This 
Is  clearly  Illustrated  In  Figure  3  where  frequency  and  cuaulatlve 
distribution  curves  are  shown  with  m  ■  7,  oT  «  0  for  oQ  -  200,  300,  400,  500, 
600. 


Once  the  Welbull  parameters  have  been  evaluated  from  a  given  set  of  fracture 
data,  the  mean  value  and  coefficient  of  variation  can  be  calculated  from  the 
respective  formulae  (Kef  4) 


a  «  Oj  +  o0  T  (1  +  1/m)  , 

C  -  ✓  r  (1  +  2/m)  -(!*(!  + 


Here  T(z)  Is  the  gamma  function  of  z.  Simplifications  occur  when  o~  »  0  and 
Equation  1  reduces  to  the  two  parameter  distribution.  The  analysis  of 
fracture  data  based  on  this  assumption  usually  yields  a  values  which  are 


greater  than  5  and  this  results  In  the  approximate  formulae 


o0  1  -  0.57  , 


C  -  1.27 

m  +  0.55 


(7) 


From  the  last  equation  It  Is  seen  that  for  a  two  paraaeter  Welbull 
distribution  the  Welbull  modulus  a  is  an  Inverse  measure  of  material 
variability. 


3  EVALUATION  t)F  THE  WEIBULL  PARAMETERS 

The  evaluation  of  the  Welbull  parameters  can  be  achieved  by  curve  fitting, 
the  method  of  moments,  or  the  maximum  likelihood  technique.  Curve  fitting 
Includes  linear  plotting  and  non-linear  least  squares  estimation  (Refs  2,5). 
For  the  two  parameter  distribution  there  Is  little  to  choose  between  the 
above  methods  and  reliable  results  can  be  achieved  In  each  case.  For  the 
three  paraaeter  distribution  the  choice  of  method  Is  important  and  111 
conditioning  can  arise  for  some  numerical  procedures.  In  this  section  a  new 
approach  Is  presented  that  Is  based  on  a  non-linear  least  squares  analysis 
which  appears  to  offer  Improved  consistency  over  previously  documented  least 
squares  techniques. 

In  Equation  1  the  Welbull  parameters  refer  to  population  values.  For  most 
engineering  applications  only  small  samples  of  data  will  be  available,  so  the 
parameters  calculated  will  be  estimates  of  the  true  parameters.  As  the 
sample  size  increases,  these  estimates  will  approach  the  population  values. 
Confidence  limits  associated  with  the  small  sample  estimates  are  clearly 
Important  and  this  topic  Is  discussed  elsewhere  by  the  authors  (Ref  1). 
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3.1  Plotting  Position 

If  the  fracture  data  are  arranged  In  ageending  order  so  that 

0  -  °1  •  °1+1  >  °1  »  1  *  1»2 . .  •  <8> 

an  estimate  of  the  probability  value  aasoclated  with  each  nay  be 
calculated.  The  true  rank  Is  unknown  but  Is  distributed  for  the  1th 
observation  according  to  the  density  function  (Ref  6) 

g  (Pi>  -  n  i  (Pt)1_1  (1  -  Pi)""1  •  (9) 

(1  -  1)  I  (n  -  1)  ! 

In  practice,  the  mean  or  expected  value  Pi  of  the  rank  Pi  Is  often  used  where 
P1  “  £  pi  8  <pi>  d  pi  <10> 


which  simplifies  to 


-  1  .  (11) 

n  +  1 

However,  computations  have  shown  that  this  assumption  leads  to  biased 
estimates  of  the  Welbull  parameters  (Ref  7)  resulting  from  Inaccurate 
representations  of  the  fracture  data  In  the  lower  probability  region.  An 
alternative  approach  Is  to  choose  the  median  of  the  distribution.  These  are 
calculated  by  solving  the  Integral  equation 


8  <pi)  d  pi 


0.5 


(12) 


to  find  a  value  of  the  median  rank,  P^  so  that  the  associated  cumulative 
distribution  function  Is  0.5.  For  n  less  than  20  the  Pj  value  is  obtained  by 
solving  Equation  12  numerically.  Results  for  sample  sizes  of  n  ranging  from 
1  to  20  are  tabulated  In  Reference  6  which  also  gives  the  approximate 
formula 


1  -  (1  -  In  2)  -  (2  In  2  -  1) 

- U - 


(13) 


for  use  when  n>20.  Adopting  these  values  for  P^  leads  to  a  better 
representation  of  the  data. 

3.2  Calculation  of  the  Welbull  Parameters 


Once  the  median  ranks  are  known  the  Welbull  parameters  may  be  calculated. 

Many  techniques  have  been  presented  for  evaluating  these  quantities  end  some 
of  the  methods  are  outlined  In  a  paper  by  Schneider  and  Pilazotto  (Ref  5). 
Most  methods  fall  Into  two  categories;  a  trial  and  error  graphical  approach 
In  which  the  data  are  linearised  or  a  least  squares  minimisation  based  on  the 
probability  curve.  The  former  method  normally  uses  a  minimisation  scheme 
based  on  the  differences 
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Si  -  Xnln  [1/(1  -  P(a1)]  -  lnln  [l/(l  -  Pt)]  .  (14) 

In  the  latter  method  use  Is  aade  of  the  difference  scheme 

«i  -  P(0i)  -  Pt  (15) 

where  P(0j)  Is  the  value  of  P  calculated  from  Equation  2  using  c^,  and  P^  Is 

the  selected  plotting  position.  Both  schemes  rely  on  a  minimisation  with 
respect  to  probability  values  and  computations  have  shown  that  significant 
changes  in  the  calculated  parameters  may  be  caused  by  small  changes  In  the 
data.  This  Is  particularly  noticeable  for  the  threshold  stress  which  In  some 


cases  can  have  a  calculated  value  outside  the  valid  range. 

0  <  oT  <  .  (16) 

More  reliable  calculations  can  be  made  for  the  Welbull  parameters  If  a  least 
squares  analysis  Is  adopted  which  Is  based  on  the  difference  scheme 

61  "  0  (pl>  “  Oi  “  -  ot  .  (17) 

In  equation  17  satisfies  the  transformed  Welbull  relation 

-  0T  +  °0  [In  (1/(1  -  P^)]17®  (18) 

for  the  selected  plotting  position  P^.  Minimising  the  sum  of  the  square  of 
the  differences 

A  -  £  («i)2  (19) 

1-1 


with  respect  to  each  Welbull  parameter  leads  to  the  least  squares  equations 

(20) 


n  °T  + 

n 

o„  Z  V 
°  1-1 

1/m  _ 

1 

n 

£  Oj  • 
i-1  1 

n 

c»»<  £ 

W.1/b  + 

n 

o  Z 

w  27“  -  z 

T  i-i 

1 

0  1-1 

1  1-1 

n 

Ot  £ 

T  1-1 

w/7®  In 

wi  + 

o0  !  w^7® 
0  1-1  1 

.  v  1/m 

?i  » 


1/m 


where 


Wj  -  In  (1/(1  -  Pt» 


(21) 

’(22) 

(23) 


Equations  20  and  21  say  be  solved  simultaneously  to  express  o0  and  Of  In 
terms  of  the  Welbull  modulus  m,  le 


n  I  a,  W.l7“  -  Z  W.17®  I  o, 

°o  _  1*1  —  1  _  1-1  1  _  w  1  • 

r  u .2/m  _  ^  w  1/m  £  H  1/s 


(24) 


n  £  W 
1-1 


1-1  1  i-1 
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!  o,  ?  W.2'®  -  §  W<’ 

-  1-1  x  i-1  1  1-1  x 


l/« 


TT 

E  "1 
1-1  1 


TT 

£ 

1-1 


n 

£ 

1-1 


a  U  1^n 
°i  W1 


—TT 

£ 

1-1 


2/m  _ 


1/n 


1/m 


(25) 


Substituting  for  oQ  and  oT  Into  Equation  22  leads  to  a  transcendental 

equation  for  m.  Although  Its  final  form  Is  too  complex  to  be  written  down  In 
this  text  the  resulting  equation  may  be  readily  solved  numerically  on  a  high 
speed  computer*  In  the  root  solving  routine  an  Initial  value  of  m  -  1.0001 
can  be  chosen  to  start  the  Iterative  solution  procedure.  For  values  of  m<l 
the  Welbull  equation  degenerates  Into  a  physically  unrealistic  form  for  the 
studies  under  Investigation. 

The  above  analysis  has  been  developed  on  the  assumption  that  the  threshold 
value  o<f  Is  unknown.  If  the  threshold  value  can  be  prescribed  certain 
simplifications  occur.  For  analytical  reasons  It  is  convenient  to  Introduce 
the  difference  notation 


-  In  (o  (Pt)  -  Oj)  -  In  (o1  -  oT) 


(26) 


which  via  Equation  17  reduces  to 

6^  -  In  (0j  -  Oj)  -  In  ( Oj  -  oT) 


(27) 


Performing  a  least  squares  analysis  based  on  Equation  27  gives  the  explicit 
expressions  for  the  Welbull  modulus  and  normal lsl-q  factor 


n 

l  w. 


n  j 

n  £  W .  '  ..  «* ,  ■■  / 

1-1  1  1-1  1  1-1  1 
IT 


n  i£i  Wj  In  (ot  -  oj)  ~  j£i  Wt  ^  In  <ot  -  oT) 


exp 


Z  W/  £  In  (0l 


U  N<  a* 

i-1  1  i-1 


'  Jl  ^  Jl  ^  10  (°t  ‘  °T> 


(28) 


n  I  V. 
1-1  1 


~ E - JT 

£  W.  I  W, 
1-1  1  1-1  1 


where 


Wt  -  In  (In  (1/(1  -  Pi)))  .  (29) 

4  BRITTLE  MATERIAL  CHARACTERISATION 

The  Welbull  distribution  defined  by  Equation  2  Is  applicable  to  the  analysis 
of  tensile  test  fracture  data.  However,  because  of  the  practical  differences 
previously  mentioned,  flexural  tests  on  slender  beans  are  normally  used  to 
characterise  brittle  materials. 
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When  a  a lender  bean  la  eubjected  to  bending,  the  stress  Is  uniaxial  but  non- 
uniform.  Equation  2  la  no  longer  applicable  In  Its  present  form  and  requires 
Bodlflcatlon  to  account  for  the  non-uniform  stress  field-  If  the  bean  Is 
divided  Into  N  elements,  each  assumed  to  be  subject  to  a  uniform  stress,  the 
failure  probability  of  the  assembly  Is 


P 


(1  -  P  (e)) 
1 


(30) 


»  <*> 

where  |[  (1  -  P  )  Is  the  product  of  the  Individual  reliabilities  of  the 

1-1  1 

N  elements  (Kef  8).  Applying  this  argument  to  an  elementary  volume  AV,  and 
allowing  AV  to  tend  to  zero  leads  in  the  limit  to  the  following  expression: 


p(V)  .  !  _  „p  .If 
V  V 

p(V)  -  o  L 


0>0j  , 

0<0T  . 


(31) 


A  similar  result  can  be  derived  for  area  defects 


P(A) 

p(A) 


0>Oj  , 

o<oT  . 


(32) 


Here  V  and  A  respectively  denote  the  volume  and  surface  area  of  the  beam  over 
which  Integration  is  performed.  For  simplicity  the  symbols  m,  oj  and  oQ  have 
been  used  In  both  Equations  31  and  32,  but  the  parameter  values  will  differ 
depending  on  the  particular  type  of  defect  being  Investigated. 


If  a  slender  beam  Is  subjected  to  a  constant  bending  moment  M,  the  bending 
stress  Is  given  by 

0(b)  .  HZ  (33) 

I 


where  I  Is  the  'moment  of  inertia'  of  the  beam  cross  section  about  the 
neutral  axis,  and  y  Is  the  distance  of  a  point  In  the  cross  section  from  that 
axis.  A  convenient  way  of  producing  the  above  stress  state  Is  to  subject  the 
beam  to  four  point  loading  shown  dlagrammatlcally  In  Figure  4.  For  the 
purpose  of  this  analysis  failure  loads  are  only  recorded  when  fracture  occurs 
In  the  central  span  where  the  bending  moment  Is  constant. 


Considerable  simplifications  can  be  made  to  the  theory  If  rectangular  beams 
arc  selected  for  testing.  In  this  case  the  volume  and  surface  Integrals 
appearing  In  Equations  31  and  32  can  be  evaluated  In  closed  form  to  yield 
(Ref  9) 


Ft(v>  -  0 


Oi<b>«fl* 
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In  the  above  equations  b  and  h  respectively  denote  the  breadth  and  height  of 
the  beam  and  o^*1)  Is  the  1th  maximum  bending  stress 


°i(b) 


(36) 


For  analytical  reasons  It  Is  convenient  In  Equation  34  to  define  the 
auxiliary  variables 


m  “  m  +  1  , 

1 

o0  -  (2  (m  +  1)  o( ,m)“+1  ,  (37) 

* 

Oj  *  Otp  a 


Equation  34  can  then  be  rewritten  as 

0l(b>>sT . 

(38) 

P<V>  -  0  . 

and  rearranged,  for  a  particular  bend  stress  to  give 

°l<b)  *  «T  +  *o  [°i(b)  131  (l/(1  ‘  P1<V)»]  l/m  (39) 


When  Equation  39  la  compared  with  Equation  18  it  Is  seen  that  they  are 
Identical  in  form  provided  that  the  term  W  defined  by  Equation  23  la 
replaced  by 

Wj  -  In  (1/(1  -  Pj(V)))  (40) 

* 

It  therefore  follows  that  the  solution  for  the  auxiliary  parameter  set  n. 


10 

UNCLASSIFIED 


UNCLASSIFIED 


UNCLASSIFIED 


than  unity  and  noting  that  for  all  caaea  a  la  greater  than  unity,  the  ten 
will  always  be  less  than  0.5.  Under  these  circumstances  It  Is  possible  to 
arrive  at  values  for  the  Welbull  parameters  by  making  use  of  the  Iterative 
equations 


o/W  - 


£i(k+l) 


Vk+1)  ♦  *0 


h 

b  (m<k>  +  1) 


1  In  (1/1(1  -  Pt(A) 

1  +  El(k+1> 


*  (k+1) 
k  -  1,2. 


(«) 


Comparison  of  Equation  4S  with  Equation  18  shows  that  they  are  of  the  seme 
form.  Hence  the  auxiliary  parameter  set  may  be  obtained  using  Equation  45 
with  Equations  20  to  25  when  In  those  equations  Is  replaced  by 

H1  “  X  In  (1/(1  -  P^*)))  .  (46) 

1  +  et(k+1) 

*  *  * 

The  Initial  values  m(l),a0(l),o.j(l)  required  to  start  the  Iteration  are 
obtained  by  solving  Equation  43  with  set  equal  to  sera.  The  Iterative 

process  la  tenlnated  when  the  relative  difference  between  each  parameter  Is 
less  than  some  prescribed  tolerance.  Computations  have  shown  that  this  is 
usually  achieved  In  less  than  ten  Iterations  for  a  prescribed  tolerance  of 
1.0x10“ 5 .  The  parameters  of  the  uniaxial  distribution  are  then  obtained  using 

the  final  values  of  m(k),  *0^k)  and  In  Equation  42. 


4.1  Size  Effect  Relation 


The  average  tensile  failure  strength  "a  associated  with  volume  or  surface 
flaws  depends  on  the  volume  or  surface  area  of  the  component  under  test. 

This  phenomenon  Is  known  as  the  'size  effect'.  For  the  flexural  tests 
described  previously,  the  volume  and  surface  area  under  consideration 
corresponds  to  the  region  between  the  two  inner  knife  edges.  It  can  be  shown 
that  the  'size  effect'  Is  governed  by  the  relation  (Ref  8) 


where  Oy^  and  Oy^  are  the  respective  strengths  of  volumes  Vj  and  V2.  Using 

this  relationship  a  unit  volume  average  tensile  failure  strength  ov  can  be 
introduced  and  related  to  any  other  volume  strength  ov  by  the  equation 
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°v  ~  “T  -  /V\1/B  .  (48) 

°V  _  °T  '  I 


Although  v  has  the  numerical  value  of  unity  It  la  retained  ao  that  the  above 
equation  has  the  correct  dimensional  form.  Eliminating  of  from  Equation  2  In 
favour  of  oy  using  Equations  4  and  48  results  In  the  modified  Wslbull 
equation 


P  -  1  -  exp 

P  -  0 


r{(l  +  l/m)}“  V  0  ~  °T  \a 
V  \°y  ~  °T/ 


oxjj  , 


0<0T 


(49) 


In  this  expression  the  failure  probability  of  a  component  of  volume  V  which 
is  subjected  to  a  uniform  uniaxial  stress  o,  Is  given  In  terms  of  the 
characteristic  quantities  m,  oT  and  oy.  A  similar  relationship  may  be 
derived  for  surface  area  defects. 


5  MULTIAXIAL  “FAILURE  THEORY 


When  a  structure  Is  loaded,  the  state  of  stress  at  a  point  will  In  general  be 
multlaxlal  and  can  be  uniquely  defined  In  terms  of  the  three  principal 
stresses  olf  o2,  o3.  For  a  given  set  of  loading  conditions  these  stresses 
vary  In  magnitude  and  direction  as  a  function  of  position  within  the 
structure.  By  treating  the  loaded  component  as  an  assembly  of  elements.  It 
can  be  readily  shown  by  considering  the  reliability  of  each  element  that  the 
extension  of  Equation  49  to  a  non-uniform  uniaxial  stress  field  Is 


P  «  1  -  exp 


P  -  0 


{r(i  +  i/»>}"  v  ;v  1°  ~  °T 


°v  -  °T 


dV 


OXJj 


(50) 


o<o, 


T  * 


It  Is  not  possible  to  extend  the  above  equation  rigorously  to  Include  a  non- 
uniform  multlaxlal  stress  state;  Instead  It  Is  necessary  to  postulate  a 
criterion  of  failure  which  will  describe  this  condition.  Many  different 
theories  have  been  presented  for  various  materials  (Refs  10,  11),  and  It  is 
the  author's  view  that  It  la  unlikely  that  a  single  theory  will  emerge  that 
will  describe  the  behaviour  of  all  materials.  Individual  theories  will  need 
to  be  developed  and  assessed  for  the  particular  materials  under 
Investigation.  On  the  basis  of  this  philosophy,  a  theory  la  presented  here 
which  gives  reasonable  results  for  a  range  of  graphites  used  In  rocket  nozzle 
applications  in  the  UK. 


As  a  working  model  It  la  assumed  that  the  reliability  of  an  elemental  volume 
la  calculated  by  considering  tha  effect  of  each  of  the  principal  atrasaes 
Independently.  This  la  the  criterion  of  independent  action,  originally 
proposed  by  Stanley  at  al  (Raf  8).  With  this  assumption  the  reliability  of 
each  elementary  volume  due  to  each  principal  stress  in  turn,  and  hanca  that 
of  the  entire  structure,  can  be  assessed.  A  generalisation  of  Equation  50, 
which  la  based  on  this  assumption  and  Is  applicable  to  anisotropic  materials 
subject  to  both  compressive  and  tensile  stresses,  is  given  by 
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p  ■  i  -  exp  -  {r<i  +  i/*}*  v  .  sy 
P  -  0  ^ 


°i>°T 

°1<0T 


(T) 

(T) 


°1>0T 

V0T 


(C) 

(C) 


(51) 


where 

Sy  -  l  /  I  \  “  dV  , 

1-1  V\V^  "  ar1(0i)j  T" 

Ot^Oi)  -  oTj(t>  ,  aVi(0l)  -  oVi<T>  for 

“r^i)  ■  °t1(C)  •  %(°i>  "  5Vl(C)  for  °i>ar1(C)  * 


(52) 


In  Equation  SI  the  summation  i»l,2,3  refers  to  the  principal  directions-  The 
quantities  oT  and  Oj  1-1, 2, 3  are  the  threshold  strengths  In  the 

three  principal  directions  for  tension  and  compression  respectively. 
Similarly,  the  unit  volume  uniaxial  strengths  in  the  principal  directions  for 
tension  and  compression  are  oy^  '  and  °vj  »  The  compressive 

stresses  and  strengths  In  these  equations  are  negative  quantities. 


To  utilise  Equation  51  the  values  of  m,  0  and  need  to  be  known  In  the 
direction  of  the  principal  stresses  at  all  points  throughout  the  material 
volume.  Determination  of  these  quantities  for  all  possible  orientations 
would  require  extensive  mechanical  testing.  It  Is  therefore  Important  from  a 
practical  point  of  view  to  Introduce  some  simplifying  assumptions  In  order  to 
reduce  the  number  of  strength  tests  needed.  For  a  particular  grade  of 
graphite  It  has  been  shown  (Ref  12)  (by  a  Wei bull  analysis  of  flexure  data) 
that  the  'Teibull  modulus  Is  independent  of  orientation.  It  has  also  been 
found  that  the  value  of  the  unit  volume  average  strength  oy  for  an 
orientation  defined  by  the  direction  cosines  (lj,  12,  I3)  is  given  by 


where  oy  ,  1-1, 2, 3  are  the  respective  values  of  ov  in  the  principal  material 

directions.  It  Is  postulated  that  when  a  three  parameter  Welbull 
distribution  Is  utilised  (with  o^X)),  Equation  53  remains  valid.  Past 
experience  has  shown  that  the  average  strength  for  a  set  of  data  Is 
approximately  constant  whether  it  Is  calculated  from  the  results  of  a  two  or 
three  parameter  analysis.  It  may  be  possible  to  assume  that  for  the  three 
parameter  case  the  Velbull  modulus  Is  Independent  of  orientation,  and  that  0* 
Is  described  by  a  strength  ellipsoid  similar  to  that  of  "oy.  The  variation  of 
0-  la  only  an  hypothesis  and  its  exact  variation  may  be  difficult  to  define 
alnce  Its  value  Is  determined  effectively  by  extrapolation  which  can  be  very 
sensitive  to  small  changes  In  the  date. 
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Equations  SO  to  S3  have  been  developed  assuming  volume  defects.  The 
corresponding  expressions  for  surface  defects  are  obtained  by  reformulating 
the  analysis  In  terms  of  the  surface  integrals. 

6  CONCLUSION 


This  report  describes  a  new  technique  for  evaluating  the  parameters  of  a 
Welbull  distribution  from  a  set  of  uniaxial  fracture  data.  An  extension  to 
this  technique  allows  characterisation  of  a  brittle  material  from  flexural 
tests  on  beams  In  which  the  stress  field  Is  uniaxial  but  non-uniform.  Using 
this  method,  the  characteristic  Welbull  parameters  for  either  surface  or 
volume  defects  can  be  obtained.  In  order  that  multlaxlal  stress  states  can 
be  analysed,  the  theory  of  Independent  action  Is  then  Invoked.  This  failure 
theory  makes  use  of  either  surface  or  volume  defect  Welbull  parameters  In  the 
failure  probability  calculation.  The  method  of  characterisation  and  failure 
prediction  using  the  theory  of  independent  action  are  verified  experimentally 
In  another  report  by  the  authors  (ref  1).  In  Reference  1,  experimental 
results  are  analysed  and  methods  are  developed  which  enable  the  probabilistic 
design  technique  to  be  applied. 
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FIG.  1 


FIG.  1  WEIBULL  DISTRIBUTION  FUNCTIONS  AND  CUMULATIVE 
PROBABILITY  CURVES  FOR  m  =  3,  go  =  1  AND  VARYING 
gr  VALUES 
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FIG.  2 


IG.  2  WEI  BULL  DISTRIBUTION  FUNCTIONS  AND  CUMULATIVE 
PROBABILITY  CURVES  FOR  ^  =  0.  =  1  AND  VARYING 

m  VALUES 
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FIG.  3 


FIG.  3  WEIBULL  DISTRIBUTION  FUNCTIONS  AND  CUMULATIVE 
PROBABILITY  CURVES  FOR  m  =  7,  gT  =  0  AND  VARYING 
qo  VALUES 
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FIG.  4 


W/2  W/2 


FIG.  4  DIAGRAM  SHOWING  FOUR  POINT  FLEXURE  LOADING 
AND  THE  RESULTING  BENDING  MOMENT  DIAGRAM 


UNCLASSIFIED 


REPORT  DOCUMENTATION  PACE 


(Notts  on  co^lction  ovsrlsaf) 

Overall  security  classification  of  shoot  . .  . 

(As  for  os  possibio  this  shoot  should  contoin  only  uncloss  if  led  information.  If  it  is  necessary  to  oator 
elossifiod  information,  the  bos  concerned  must  bo  marked  to  indicate  the  classification  eg  (f),  (C),  or  (S)). 


1.  DftIC  Reference  (if  known) 


2.  Originator’s  Reference 

Report  16/84 


3 .  Agency  Reference 


A.  Report  Security 
Classificetion 

UNCLASSIFIED 


3.  Originator's  Code 
(if  known) 

7699000H 


6.  Originator  (Corporate  Author)  Name  and  Location 

Ministry  of  Defence 

Royal  Armament  Research  and  Development  Establishment 


5A.  Sponsor ing  Agency's 
(Code  (if  known) 


6A.  Sponsoring  Agency  (Contract  Authority)  Name  and  Location 


7.  Title  Application  of  Modified  Three  Parameter  Weibull  Distributions  to  Brittle 
Material  Design 


7A.  Title  in  Foreign  Language  (in  the  case  of  translation) 


7B.  Presented  at  (for  conference  papers).  Title,  place  and  date  of  conference 


B.  Author  1,  Surname,  initials 

Mar get son,  J. 


9A  Author  2 

Cooper,  N.R 


9B  Authors  3,  6. 


10.  Date  pp  ref 

February  1986  20  12 


111.  Contract  raber  and  Period 


12.  Project 


134  14  Other  References. 


15.  Diatribution  statement 


Descriptors  (or  keywords) 

I^^Brittle  fracturing  Mathematical  prediction 

^  Brittleness  '  Weibull  density  functions 

Design  criteria  Failure 


Axial  stress 


Abstract"^  The  three  parameter  Weibull  distribution  is  used  as  a  basis  for 

probabilistic  design  with  brittle  materials.  Beam  fracture  data  are  analysed  by 
a  new  technique  to  obtain  the  characteristic  Weibull  parameters  which  are  then 
used  in  conjunction  with  the  theory  of  independent  action  to  predict  failure  in 
a  multiaxial  stress  state.  ►  ,  ,,„r  - 


